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Abstract. We study a finite Larmor radius model used to describe the ions distributions in the core of a toka- 
mak plasma, that consist in a gyro-kinetic transport equation, coupled with an electro-neutrality equation. Since 
the last equation do not provide enough regularity on the electric potential, we introduce a simple linear collision 
operator adapted to the finite Larmor radius approximation. Next we study the two-dimensional dynamics in 
the direction perpendicular to the magnetic field and prove thanks to the smoothing effects of the collisions and 
of the gyro-average the global existence of solutions, as well as short time uniqueness and stability. 
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1 Introduction. 

The model studied in that article describes the density of ions in the core of a tokamak plasma. In such highly 
magnetized plasma, the charged particles have a very fast motion of gyration around the magnetic lines, called 
the Larmor gyration. A good approximation is then to consider that the particles are uniformly distributed 
on gyro-circle, parametrized by their gyro-center, and Larmor radius r^ (that is proportionnal to the speed of 
rotation u, and in our article, we will forget the physical constant and write r^ = u). The models obtained 
in that new variables are kinetic in the direction parallel to the magnetic field lines, and fluids (precisely a 
superposition of fluid models) in the perpendiciilai' direction. For rigorous dcrivatign of such models and more 
complete discussion on its validity, we refer to |tJ and our previous work [^, m which the derivation is perform 
from a Vlasov equation in the limit of large magnetic field. 

Such gyro-kinetic models are usually closed by an electro-neutrality equation, that as usual provide very few 
regularity for the eletric field, so that the well-posedness of gyro-kinetic models is, at least at our knowledge 
unknown. In this article, we add a "gyro-averaged" collision operator to the model and study the dynamics in 
the directions perpendicular to the field only. 

Let us now describe our precise model. The ion distribution function J{t,x,u) in gyro-coordinates depends 
on the time i, the gyro-center position cc € T^ and the velocity of the fast Larmor rotation u E M+ (which is also 
proportional to the Larmor radius). The electric potential $ depends only on {t,x). They satisfy the following 
system of equation on fi = T^ x R+ 

% + ( J°V,$)^ • V,/ = Puduf + 2/3/ + ^ (a,/ + -duiuduf) 
ot \ u 

($-$*, HT){t, x)^T (pit, x) - 1) 
p{t,x)^ {fl,f{t,x,w)2TTwdw) 



f{o,x,v) ^ fi{x,v) y{x,u)en 

where /3 and v are two positive constant, p is the density in physical space, T is the ion temperature. 
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We also used the notation 6^ = (— &2, bi), for any vector b = (61, 62) of M^. 

|Gh eHauNou09 

That model without the Fokker-Planck operator {v = 13 = 0) was studied in ourprevious work p - to 
which we refer for an heuristical derivation pf„the electro-neutrality equation (|IA^Il - and is used by physicists for 






simulation, by instance in Gysela code \5 . Here we just mention that ^li.z]) is obtained in a close to equilibrium 
setting, with an adiabatic hypothesis on the distribution of the electrons Ue = n^e—^ ~ no ( 1 + ^ ) i ^-nd an 
hypothesis of adiabatic response of the ions on the gyro-circle which gives rise to the $ * Ht term. As usual in 
quasi-neutral equation, we have no good a priori estimates on the regularity of _E = — V<i>. 

lag ■ iry^"'^P'^" 

Remark that even if that equation (|ll.l|) is derived frcim aYlasov model (A rigourous derivation of a more 
general iD model is performed for fixed field E in section |Ef, it is of "fiuid" nature. In fact there is no transport 
in the remaining of the velocity variable w, and the position of the gyro-center is transported by the eletric drift 
(J° £')-'-. So that the equation is similar to the 2D Navier-Stokes equation written in vorticity. More precisely, 
we have a family of fluid model depending on a parameter u, which are coupled thanks to diffusion in the u 
variable, and by the closure used for E described below. 

Moreover, we will prove in the following that thanks to gyro-average J^, equation has the same regularity 
than the NS2D equation in vorticity. In fact, the force field J^V x^ belongs naturally to H^ ii f G L^ with some 
weight. That is why we obtain the same result that are known about the NS2D equation : global existence and 
short time uniqueness and stability. However, our model present an additional difficulties which is the lose of 
regularity for small u. In fact, for small value of u the H^ bound (in x only) of J,°Va;<I> explodes. 

To state our reuslts properly, we will need the following definitions and notations : 

• In the sequel, the letter C will design a numerical constant, that may change form line to line. Unless it is 
mentioned, such constants are independent of everything. 

• i^(^) = L'^{fl, udxdu) is the space of square integrable functions with respect to the measure udxdu. 

• We shall use various norm on T^ or on fJ. To avoid confusion, we will use the following convention. All the 
norm performed on the whole fi will have their weight with respect to u as additional indice. By instance 
II • II 2;™, II ■ llffi ■ AH the norms without any indices are norm on T^ only. 



For any weight function k 



H^ 



the norm 



is defined for any function / on fi by 



l/ll 



2,fe 



\F{-,u)\\2k{u)du 



The most usefuU weights will be m{u) — 27ru(l -t- u^) and ■m{u) = I + u^. 

We change a litlle bit the duality used to define distributions in the following definition 

Definition 1.1. Using distributions with the weight u means that duality is performed as 

{fi9)u = / fgdxgdv udu. 



This definition may seem a little artificial because the simple definition of derivative with respect to u, is 
not valid. Instead, 

{duf,9)u = -{f,dug)u- {-,9)u- 

u 

However, this weight respects the underlying physics {u is in fact the ID norm of a 2D velocity variable) 

and has many advantages. For instance the operator {l/u)du{udu) is self-adjoint with this weight. 



Our precise result are the following. We prove global existence under the hypothesis ||/i||2,m < +00. 



Theorem 1.1. Let J, ^'^>MsfMJ^fj}\2 m < +00. Then there exists at least one weak solution f ^ L 



oo('B + 



L'^{W^ , H^{fl)) to (|ITTl|H|li.L^|) with initial condition fi, which also satisfies for any i > 



,Lim)n 



II/WIIL + W ll(v.,a„)/||Lrfs<||/.||2,„, 
"'0 

and all the a priori estimates of the previous section (Lemmalcl. bi IIJ. bi |I5'. 'I\i if their initial hypothesis are satisfied. 



And we prove short time uniqueness and stability under the additional hypothesis ||V2:/||2,m < +00. 
Theorem 1.2. Let fi satisfy 



Il/i|l2,m + |lV:r/||2,m < +00 . 

n fim : p|rad_xf Jflg : FY_ - 

'hen thenositive time r* defined in Leninia p.7\ is such that the weak solution to (|ll.l|) -|l 
Tl IS unique on [0, r*] . 
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Moreover, that solution is stable jm that interval of time in the following sense. Assume that (f")rif=N is a 
' Jftini - existence ■' . . y'lti-a 

family of solutions given by theorem \\l . 1\ with initial conditions /" satisfying 



lim ||/r-/dl2,™ = 0, 



sup 1 1 /r 

nSN 



\li^(l^) < +00 , 



Then 



lim sup 



irw-/(i)ii2,™ = o. 



l2,m ■ 



This local result has some more consequence when relating it to the bound on zy/„ ||Va:/||2 m^t < \\fi\\t 
C(T)||/i||2 satisfied by any solution in the sense of ill. 1 1 Ihe last bound implies that ||Va;/||2,m is almost surely 
finite. The local result implies more : that the norm of gradient may blow up only on a closed and negligeable 
set, of 4/5-capicity zero 

lag : gymFP^n 

In the next section the diffusive operator of (|ll.ll) is rigourouslv derived from a linear Vlasov-Fokker-Planck 
equation in the limit of large magnetic field. In the third section, some useful lemmas are established, proving 
regularizing properties of the gyro-average, global preservation of 



preservation of the u(l 



ic weighted norm of /, the short time 
■ M'j-moment of Vxf by the system (pi.i|) - ipi.:!il) ). and, controlling the electric potential by 
the physical density. This allows to prove the global existence (Theorem 111. l|l of solutions to the Cauchy problem 
in the fourth section and their short time uniqueness and stability (Theorem u.Z\] in the fifth section. Finally 
some useful properties of the first Bessel function j" are proven in the appendix. 
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2 Derivation of the gyro-Fokker-Planck operator 

In that section, we rigorously justify the form of the Fokker-Planck appearing in the right-hand side of ((ff 



YroFP2D 



The usual collision operator for plasmas is the nonlinear Landau operator originally introduced by Landau \T~ 
Because of its complexity, simplified collision operators liavc been introduced. An important physical litterature 
exists on the subject, also in the gyro-kinetic case (See \^ and the references therein). In this paper we choose 
the simplest possible operator possible, namely a linear Fokker-Planck operator. The reasons of this choice are : 

- Its simplicity will allow to focus on the other difficulties of the model, 

- The fact that physicists studying gyro-kinetic models for the core of the plasma mainly assume that the 
dynamics stays close to equilibrium, in which case a linear approximation of the collision operator is relevant. 

- The aim of the paper is not a precise description of collisions. In fact, even if they exists in tokamaks, being 
needed to produce energy, their effect is small compared to the turbulent transport. However, we are interested 
by their regularizing effect, since the electro-neutrality equation (|ll.:^|) do not provide enough regularity to get a 
well-posed problem. This is a major difference to the Poisson equation setting. 

We start from a simple model for a 31? plasma, i.e. a linear Vlasov-Fokker-Planck equation with an external 
electric field, an external uniform magnetic field and linear collision ^d drift terms, and obtain in the limit 
of large magnetic field a 3D (in position) equation analog to (|ll.l|) . in particular, we show that a usual linear 
Fokker-Planck term on the speed variables turns into an equation with diffusion terms both in space and Larmor 
radius variables in the limit. 

Precisely, for any small parameter e > we study the distribution fe{t, x, v) of ions submitted to an exterior 
electric field E{t,x) (independent of e) and an uniform magnetic field B^ — (1/e, 0,0). We also model collisions 
(with similar particles and the others species) by a simple linear Fokker-Planck operator. To avoid any problem 
with possible boundary collisions, which are really hard to take into account in gyro-kinetic theory, we assume 
that (x,v) gT^ X M.'^, where T'^ is the 3D torus. When the scale length of all the parameters are well chosen (in 
particular the length scale in the direction perpendicular to the magnetic, field ^hould be chosen of order e times 
the length scale in the parallel direction, we refer to our previous work g for more details on the scaling), the 
Vlasov equation /^ satisfies is 



ridau 
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-J-+v^^d,^^f + E-VJ + -iv^-VxJ- 



V„^/) = div„(/3t;/,) + i.A„/,, 



(2.1) eq:vlares 



where /3, v are two positive parameters, the subscript || (resp. _L) denotes the projection on the direction parallel 
(resp. on the plane perpendicular) to -B, and the superscript _L denotes the projection on the plane perpendicular 
to B composed with the rotation of angle 7r/2. In others words li v — {vi,V2,v^), 



V±_ = (ui,W2,0), 



(0,0, Vg), «^-(-f2,t'l,0). 
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The next results require the additional notation, 



27r 



27r 



Jugi^giPL^^)^ TT gixg+ue^^^Pr^e'^"^" ^^'> + v\\e\\) dfc , 



(2.2) 



which is a gyro-average performed in phase space, that will be used as an initial layer to adapt the initial condition 
to the fast Larmor gyration. 



ylares 



Theorem 2.1. Let E £ L^{L ) and f^ he a family of solutions to equation (p.i|) witti initial condition fi G L 
satisfying sup^ ||/e(i)||2 < II /ill 2- Then the family /j defined by 

Mt,Xg,v) = f{t,Xg + v^,v) 

admits a subsequence that converges in the sense of distributions towards a function f depending only on 
{t,Xg^u — \v\^v..) and solution to 



dtf + v« d,J+ 4E« d,J+iJuE)^ ■ V,,/ 



/3(i;„ a,„ / + ud^f + 3/) + z. A,^J+ -duiudaf) 



O/J-O^ 



in the sense of distributions with the weight u, with the initial condition J^{f' 

Remark 2.1. The reason for the change of variables is that the 1/e-term in equation (|L^.i|l induces a very fast 
rotation in the perpendicular direction both in the x and v variables, 

v{t) = v"e''/' , x{t) =x° + «°^ + ^0e'(*/e-V2) . 

But in the gyro- coordinates this fast rotation is simply described by a rotation in v, 



vit) 



- „o^**/e 



■^gy'') ^ -^g 



Remark 2.2. The final diffusion appears in all dimensions except the Xg 



It does not 



that there is 



no regularization in that direction. Indeed, the models have diffusion in w , which after some tirrie regularize in 
the Xg direction. This mechanism is well known for the Fokker-Planck equation (see for instance [T/y -However, 
we are not able to prove this phenomena in the non-linear setting because the electric field of the model lacks 
regularity. This is the reason why we will only study the 2D model. 

Jthni ■ FPp-v-rn pheHauNou09 

Proof of Theorem llj^.i! We proved in a previous work p^ that, provided / G L and EL£z Li ( W^' ) , a subsequence 
of /e solutions of (|L!.l|l without the collision term converges towards a solution of ^l.'6\ without the collision term. 
In order to simplify the presentation, we will neglect the electric field and the parallel translation terms. To 
obtain the result in full generality, the only thing to do is to add the argument given in our previous work to the 
one given below. For the same reason, we shall also not treat the problem of initial conditions. 

So consider the above Vlasov Fokker-Planck equation without electric force field and parallel translation. 



9t/+-(f±-V,^/ + z;^-V 



,^f)^div,{Pvf)+uA,f. 
Since 



(2.4) 



,,-L 



The first step is to use the change of variables (x, w) -> [xg ~ x -\-v^ ,v 

V„-(«/) = vVj+if-vVU, 



(2.3) eq:FPlim 



eq: vlaFP 



equation (|L!.4|) becomes 



dtfe + ^v^ ■ Vvfe ^~l3{y- Vj, + 3/e - « • V^J,) + i/(aJ, + A, J, - 2V, • V,Je) . (2.5) 

By hypothesis fe is bounded in L'^{L'i ). Therefore, at least a subsequence of (/e) converges weakly to some 
J € Lf{L'^). Passing to the limit in (p 4|) . it holds that 

since all the other terms are bounded. For v — {ue^"^ , v ) where Lp is the gyro-phase, the previous equality means 
that / is indepeiiclent of the gyro-phase (in the sense of distribution and thus as a L? function) . 
Equation (p.4|) tested against a smooth function g independent of the gyro-phase writes 

j h(dtg-p{v-V,g~v-VJ^^g)-u{/\,g + ^^^^g-2Vi^-V,g)^ dxgdv^Q. (2.6) |eq:FP2 

We may also pass to the limit when e tends to zero in this equation and obtain that the same equality holds for 
/e replaced by /, considered as a function defined on T^ x R^. 
For the change of variable v = (ue*"^, w,, ), 

V„.g = (e^^dug + ie"^d^g, dy^^ ). 

Hence, for any function g independent on the gyrophase tp, it holds that 

'^v„g = dl g + -du{udug) , 

II u 

i^i ■ V.,,).g = V^^ ■ {e^d^g) = e^^ ■ Vi^dug, 
w- V„g = v^^dy^^g + udug. 

teq ; FP? 

The other terms appearing in (|b.bl) remain unchanged. Then, 



f(dtg - (3{v dy g + udug - ue'f ■ V^ 5) - i^{dl g + -5„(u9„g) 

V " II \l u 

+ Arc^^g- 2e"^ ■ Vi^dugyj dxgdv^^27rududip ^ 0. 



(2.7) 



Since / is independent of tp, performing the integration in ip first makes the term containing ip vanish. So the 
function / of the five variables {xg,u, v^^) satisfies 

fidtg - Piv^^dy^^g + udug) - i^idy g + -du{udug) + A^^^gU dxgdv^^udu = 0. (2.8) 

It exactly means that / satisfies the equation 

dtf = f3{v^^dy^^ f + udj+3,f) + V [dl f + A,^J+ -duiuduf)) , (2.9) 

I II V II ^ u ) 

If.q. FPIiin 

in the sense of distributions with weight u. It is the equation (|b!.a|) witliout parallel transport nor electric field. D 
If we look at solutions of this equation invariant by translation in the„direction of B, we exactly get the 

~ Tf.q: FPliin ' ■^ ^ 

2D-model announced in the introduction. In fact, if / is a solution of %l:6\ . tlien 

f{t,x,u) = / /(i,x,M,W||)dW|| 

i]). Such an assumption on / is reinforced by experiments and numerical simulations, where it 
is observed that the distribution of ions is quite homogeneous in a;,, . 



ct : apriori 



lem : 3Deq4D 



3 Some useful lemmas 

We prove here some a priori estimates useful for the proof of our theorem. In order to simplify the proof of some 
of the following Lemmas, we sometimes uses the following formulation of (|ll.lp with the genuine two-dimensional 
velocity variable. Denote by /(t, x, u) = /(t, x, |u|), u G M^. It is solution (in the sense of distribution with usual 
duality) of the following equation with AD in space and velocity variables 



dJ-VUJl'ui^) ■ ^J = i^(A./ + Aa/) + /3(2/ + u- Vs/). 



(3.1) 

e JFP4D teg : t ryrriFP^r) 

i]) IS a solution of (|ll.l|l . We can state for instance a precise 

Lemma in the case where 4> is fixed and smooth. 

Ipq ■ gYT-nFP9n 

Lemma 3.1. For a fixed smooth potentmL^, J is the unique solution of (|li.i|l untti initial condition fi if and 
only if f is the unique solution of lp.\\) unth intial condition fi . 

3Deq4D 
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Proof of the Lemma ll3.ll : iHe proof relies on the uniqueness of the solution to 
conservation of the radial symmetry of the solution. 



,VFP4D ILa dyzen 

(See \6 and the 



eq:VFP4D 



lem:phireg 



lem:rhoreg 



3.1 Regularizing properties of the gyro-average. 

In this section, some regularizing property of the gyro-average operato are proven. They are based on the fact 
that ,P ^ k^2 for large k (the precise bound are proved inllA]r. which implies that J^ maps i/* onto i^+a. It 
is important since the formula (|ll.:^|) giving the gyro-averaged potential in terms of the distribution / involves 
two gyro-averages, and thus a gain of one derivative for the gyro-averaged potential w.r.t. /. However, the 
regularizing properties of J^ are bad for small u, which raises difficulties. 

The first lemma of this section gives the regularity of the gyro-averaged potential in term of the potential 
<&. The second one gives the regularity of the density p in terms of the distribution /. We will need the two 
following definitions before stating it. 

Definition 3.2. Let f be a measurable function defined on 17. Denote by 

the norm with the weight m{u) = 2nu{l + u^). 

For any U > 0, let F be a measurable function defined on flij = T^ x [0, U]. Denote by 



\F\ 



Hh 




V^f\^ + \duf\^)27Tudu 



The lemmas stating the regularity of $ and p are the following. 
Lemma 3.2. For any s G M, u > and $ with 0-mean, it holds that 

^) ||J°$||ff= < mHs, 



IJ$I 



'H" 



lit) \\duJ^<^\\H^ < 



2~ 



-^11*11 s 



As a consequence, for any C/ > 0, 



\fMH}<^^um 



H2 



defrho 



Lemma 3.3. For any s > 0, if J f2'Kudxdu = 1 and p is defined by (|ll.5|) . then 

||p-l||^.+ i < 2in\\ fU^^^H^. 



(3.2) 



J laa^: phirflmn^: rhoreg 

Proof of Lemma lb. :!il and |I3!3I 
Denote by <t{k) the k Fourier coefficient of $. Then 



ij^wi^^ = E \Juikrmkr < E i*wi' = ii'^ii^^ ' 



fc=i fc=i 

"nil n /^^Tn ■ boundJ 

using the bound || J ||oo < 1 proved in Lemma [A. 1 1 It is~the inequahty i). For the second inequahty, remark that 



and use it together with ii) of Lemma lA.lI m 



' + '"' - ' '^'^<^, ker, (3.3) p^;^ 



w^l/cp 



•^u^ii' +^ = Ei^«$(fc,")P(i + ifcpr+^ 



fe/O 






fcT^O V I I 



< ^ll*li„. 



2 



n :p>iiT-R p ; 

For the third estimate of Lemma \ii.'Z[ remark that 



(a„jl]$) (fc) - a„ (j"(|A:|u)$(/c)) = |fc|$(fc)j"'(|fc|u) 



n t^^Tn ■bound J 

and use the bound Hi) of Lemma llA.il to get 



From this, we obtain 



\{du4mk)\ < \l^-^\m\ 



|5„(J°$)|U. < -^ii*ii^=+i. 



\nH}, 



The point iv) uses the previous inequahties. First remark that the norm || ||/ji is also equal to 

1 

{\\d^Fi;u)\\l. + \\Fi;u)\\l,y7Tudu\ . 

Using this formulation and ii)- Hi) leads to 

/ (||a„J°$||i2 + ||J°$||l^i)27rudw 

< 27r||$||2i/ t^LX±udu<im\m 
"'- Jo u 



Kn-Hi, = 



which gives the desired result and ends the proof of Lemma lb .1^ 

[Ltrhoreg 



1 + V2 
—I 

:phireg 



'H2 



Proof of Lemma IK. 51 
Denote by p{k) the fc-th Fourier term of p with respect to the space variable, i.e. 

p{k) = 2n f J°{\k\w)f{k,w)wdw. 
7 



By IB 



iboundJ 



\m\ < 27r 



\f\{k,w)w 
(l + u;2|fc|2)i/4 



dw. 



ale 



It follows from (|b^ that for k ^ 0, 



(l + |fcn^^|/5(fc)| < 2^^ / |/|(fc,u;)(l + |fcnfV^du; 



oo \ 1/2 / /.oo 

2/i„ „,,N/'i I I l.|2\s„,vi I „,,2\ 



< 247r / |/|^(fc,u;)(l + |fc|'')"w(l + u;^)dw; 



oo \ 1/2 

2*7r( / |/p(fc,u;)(l + |fcn^27ru;(l + w2)(iw 



dw 



(1 + ^^) 



1/2 



Hence, since p(0) = Lg /"(s^) dx = 1 by mass conservation 

IIP-1 



_H-=+2 



< 237r 



\ 



^ ( /" |/(fc,w)|2(l + |fc|2)f27rw(l + w2)dw 



< 237r I / ||/(w)|||j,27rw(l + u;^)rfu; 



1/2 



J Xsjil: rhoreg 

and Lemma lb. al ls proved. 
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3.2 Control of the potential by the density. 

Denote by Lt the operator that maps any function $ on T^ with zero mean to ^ (<& — <& *a; Ht) and by H^iT'^) 
the space of H" functions with zero mean. This section is devoted to a proo^jOtthe bpundedness of L^^ from 
iIg(T'') onto H^i^''-). Recall that in a Fourier setting (See the Appendix of [^ for more details), the operator 
Ht = I — TLt is the multiplication by 

9 /'+00 

T Jo 
Lemma 3.4. The Fourier multipliers Hxik) satisfy, 

\l-H{k)\>^-^{l-e-i^) , VfceZ2\{(0,0)}. 
As a consequence, the operator L^ maps any Hq, s G M., into itself with norm 

\\L-'\\m<CT 



1 — e T 



(3.4) 



boundL 



Remark 3.3. Lemma \\cl.4\ shows that ||i^^|| is bounded for small T, and of order T for large T, the physical 
case of interest. 

The houndedness of the spatial domain is essential. When defined on the whole space M^ rather than on the torus, 
the operator L^^ is not bounded. Its norm explodes in the low frequency range. 



■ LLaip.: boundL „ P o"i ■ boundJ 

Proof of the Lemma lb.4l l wo bounds on J (l) are used, namely one of the bounds of Lemma IIA.il tor l~> 1 and 
the following bound given by the Taylor expansion of J'^ near for / < 1, 

0<(J°(0)'<1- J, ifO<Z<l. 



sec : moment 



Consequently, 



oo 



l^H.,1 . I /-(,-<^) .-..■/...„. ^y_.^*.„ 



/•™ / 3,2 \ 2 /-OO ^ 

< 2/ il~ —^je--' xdx + V2w e'-' dx 

O 1 /"CXD 

- 4 4u;2^ 7„ 

where w = (Ifcl-v/T)^^- Now, using the bounds 2^5 < | and 

J W J U 



/•OO - 

T/; / e~^ dx < I xe"^ dx = — 



it holds that 

l-|ffT(fc)|>^(l-e-^). 

J IailL: boundL - 

This is the first claim of lemma lB.4l Ihc tunction of w in the right-hand side of the previous inequality on \H{k)\ 
is decreasing and goes from -j at to at +cx). Consequently its minimal value are obtained for large w i.e. for 
small |fc|, namely |fe| = 1. Precisely , 



sup|7?T(fc)|>^(l 



T 
l-snp\HT(k)\ > - (l-e~T 



Since the Fourier representation of L^^ is the multiplication by T(l — H{k)) ^ we obtain that in any Hq, s e M, 



-1 



T 4 



Lt \\h- = sup- — — - < 



k^o\l-H{k)\ 1-e- 
which is the desired result. 



3.3 Propagation of L^ and L'^(L'^) norms of /. 

The two following lemmas will be useful in the sequel. 

Lemma 3.5. Assume that \\fi\\2 u ^ +00. Then, any solution of (|ll.i|l and l|ll.4|) . for regular votential cb. satisfies 

yt>Q, ||/(t)||2,n<e^*||/,|J2,«. 

Assume moreover that ||/i||2,m < +00. Then any solution f satisfies 

p2<3i „ 1 

ymlm < WMln. + (2^ + /3)^^\\mh.u • (3.5) 



eq:u-momen" 



with the convention that ^-g — — 2t if (3 — 0. 



YroFP2D 



lem:disp Lemma 3.6. Assume ||/j||l2 (^4) < +oo and f is a solution of (|ll.i|) mtli initial condition fi with a regular 
potential <p. Then f satisfies 

ll/WllL^.(L-)<e(^+2'')1/.||L?„(L^) (3.6) |eq:disp 



Remark 3.4. A more careful analysis will show that 

g2/3t _ 1 



WfimiuL-) < ii/.iii^(L4) + (2^+/5)^^— ii/.iii^_(L.) 



Jlamjdisp 
but the simple estimate of Lem,m,a \\J.t)\ will be sufficient. 



J lajL: u-moment J saJFP4D ~ 

Proof of Lemma iB. a j Multiply equation (|b.i|l by f. Using the notations 



u=\u\, g{t,u) = -\\f{t, -,11)111, 
and integrating in the x yariable leads to 

dtg-iyAsg = -\\{V,,Vu)f{t,;u)\\l + l3{Ag + u-Vs9). 



(3.7) 



(3.8) 



Multiply the preyious equation by k{u), where fc is a smooth function on M^ with compact support and 
integrate in the yelocity variable u leads to 



dt [ I g{t,u)k{u)du] + / \\{V^,Vrdf{tr.u)\\lk{u)du^ 

{i'Ask{u) + 4pk{u) - P div{k{u)u))g{t, u) du. 
By approximation, this is still true for functions k with unbounded supports. For k{u) — 1, 



dt [e-^^' J git, u) du^ = -e-2^* J \\ (V,, V„)/(i, •, 



u)\\ldu< 0. 



Coming back to the ID original quantities, it means that 

Il/(i)||2,u<e^*||/,|l2,„. 

For k{u) — rh{u), then Ak — 4 and 

4m(u) — diy(m(M)M) = 2m(u) — rh' {u)u — 2 . 

Therefore, 

/ g[t,u)rh[u) du < / g[Q,u)ifi{u) du + 2{2v + P) I / g{s,u) duds 

Or in other words 

\\fmim<\\M\l,n + 2{2^ + l3) f \\.f (Minds. 

Jo 

Using the bound of equation (|b.!j|l . we get 

ll/Wll2,™<ll/dlL„ + (2^ + /3)^^ll/dlL- 



(3.9) 



eq: calc2 



with the conycntion that '^ „ "*" = 2t ii P — 0. 



J lfijIL:disp 

Proof of Lemma lb. bl in order to simplify the presentation, we will first performed calculation without justi- 
fying eyery integration by parts and diyision. But once we obtain an a-priori result, we will explain the small 
adaptation needed to make it rigorous. First, we denote 

a{t,u)= f\f{t,x,u)\Ux^\\f\\t, i{t,u)= [\ff\Vsffdx 



(3.10) 



Multiply equation (pi.ip by c>sign( f)| f |'^ and integrating with respect to x leads to 



dta = -Uiy / f\iV,,Vs)f\^dx + iyAsa + 8Pa + l3u-Vaa. 



Hence, dividing by ^/a 



ftVa = — -= < -^v—= + — -^ + 4/3Va + P^—j=- 



eq: alpha 
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Now, we multiply by m(u), integrate with respect to u 



/arh(u)du\ < —Qv 



—j=ifi{u)duA / — ^^m{u)du + 4/3 / ^/a.fh(u) du + f3 

J a 2 ./ J a 



;^—m[u) du 

2y^ ^ ' 



With the help of some integration by parts, we get that 



rh{u) du = —2 / \/a{'m{u) + u ) du 



2V^ 

Asa ~ f ',,^ o /" ""• Vsa 
— -j^myujdu = —2 / ;= — du 



2ai 



■rh{u) du 



/adu 
Thanks to that, the previous inequality simplify in 

9t / \fam{u)du ) < —Qv / —= fhiu) du -\ — 
\. J J \/a 4 



|Vsa[ 
2a i 



■rh{u) du . 



iV^al' 



■m{u) du + 2{P + 2iy) / y/adu 



Next we can estimate |Vaa| in terms of 7. In fact by Holder inequality 

V^a = Ws(lfUx)=4ffWsfdx 

iVsaP 



< 



wfJfUx) (J f\^sf\^dx] =16a7 



So that the second term in the right hand side of the previous inequality is controlled up to a constant to the 
first one. We precisely get 



dt 



nym{u)du \ < —2v 



7 



■m{u) du + 2{I3 + 2v) I yfadu 



(3.11) eq:ineqdif: 



From which we conclude easily. 

In the previous calculation, we have not justified all the integrations by part. To make the argument rigorous, 
a possibility is to choose a smooth function ^1 from R+ into [0,] such that S,{u) = 1 if u G [0, 1] and ^(m) =0 if 
u e [1, +CX)), and define for all C/ > iu{u) = ^ (^). Remark that It/^^U < l^'U and Wi'^jU < I^U. Xhgn ^^^.,, 
we performed the previous calculation with the weight rhu — rh^ij, we obtain an inequality very similar to (p.iljl 



dt I / \/a'm{j{u)du\ < —2v / —j=rhij{u)du' 



2(/3 + 2^) + ^ 



idu 



(3.12) eq:diffapp: 



from which we get ||/(<)||l2 /^4-, < e^ "^^^ II/jIIl^ c^4-), which give the desired result letting U going to 
infinity. 

The other point not rigorously justified is the division by ^Ja that may be zero. However, since we have 
a diffusion equation, it may be proved that for t > 0, a > Q everywhere. Or we can use a family of smooth 
approximation of yf-. Or we can say that a + e satisfy (p.iu(l with a additional term that has the good sign, 
so that it will satisfy (|b.l:^|) . and we will obtain the desired inequality letting e going to zero and then U going 
to infinity. It is well justified since the maximum principle applies there so that any solution with non-negative 
initial condition remains non-negative. 

3.4 Short time estimate of the m-moment of Vxf- 



The following lemma provide is central in the proof of the stability and uniqueness of the solution for short time. 

teg : gyT tei yP t miata 

Lemma 3.7. Assume that f is a solution of the system (|ll.l|) - (pi.4l) satisfying initially \\'^xfi\\2,m < +00. Then 
there exists a constant C* and a time t* depending on {T^v, || V:E/i||2,m), such that 



|V,/,| 



(v„a„)v,/|j2 di<c* 
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J aaif lecneutr 

We also mention that the resuh is true if tlic„dcfinition of $ in (111.:!; I) is replaced by another definition which 
stiU satisfies the bound given in Lemma lb.al ancl ltj.:!il Vrecise bounds by below for t* are given at the end of the 
proof (only in the case /3 = 0). 

Proof of Lemmalia.YI: Wc take the a;-gradient of equation (|b.l|) . written in 21? in u (with u = |u|), and obtain 



dt^.J - ^t{J»^iJ = /3(2V,/ + u ■ Vs(/)) + ^A,,„-(V,/) - V,(V^( J:;<I')V,/. 

If we now multiply by ^Wxf on the left and integrate in x, the function h defined by h{t,u) — 2 J \'^xf\'^ dx 
satisfies, 



dtg{u) ^ l3{Ag{u) + u ■ V u9{u))vAr,g{u) -v\\Vx,iNxI\[ 



'VxfVx{Vi;jy)Vxfdx. 



(3.13) eq:gevol 



We may also multiply this equation by m(u) = (1 + v?) and integrate it in u. We obtain after that 
1 



75t||V./||U + t^||V,,„-V,/||2, 



(4Z. + 2/3)11 V,/||2 „o -yy 'vj{Vx{VJ;Jl^))vj7n{u)dx 



du. (3.14) 



To go on, we need to understand a little better the matrix M(t,x,u) = Va;(V;^J°$). First remember that 
$ = L^ {p— 1), and then remark that from there definition, j" and L~^ commute with deri'^tion in ic. So that 
our term may be rewritten M = J^L^^ iy ^iy x P)) ■ Using the bound of the Lemma IbV^I and IHTJI we obtain that 



Vu>0, \\M{t,u)\\H^^\\JlL-\yx{yip))\\m<'^-^\\p-l\\„^<^\\f\\L^Hi^). 



'H^ 



I pq : Pf vol 

Moreover, the iJ„j-iiorm of / appears in the right-hand side of (|b.lj|) . So that we may use it to control M. 

With a control on the H^ norm of M, we do not get an infinite bound on M, like ||M||oo < C'||/||i?2. In 
that case, we will be able to use a classical tool to conclude. But this is almost true, we are in a critical case 
{d = 2 and p — 2) for the Sobolev imbeddings, but we still know that the square of M is exponentially integrable. 
Precisely, since M is of average 0, we have for all u > the following Triidinger inequality 



"hi dx <2 



We refer to pTor a proof of that result. To estimate J J *V/MV/(1 + u^) dxdu, we first perform the integral 
in X. For this, we apply the inequality 



ab < e°- - h + hlnh , \/a,b>Q, 



\M\ 



I2 — ) , |L,^i2 )• The previous inequality comes from the Legendre transform of e°, as the 



to (a, 6)- M6||Mr^,y Mtv/lt 

Young inequality, so that our application will be a log-exp analog of the Holder inequalities. We obtain 



||AfV/||2 = 6||Af||^i||V/||2 




\M\ 



mm. 



|v/T 
|v/||i 



dx 



< 



< 



;^ll/lk,(..llV/|b 2-1 + y^ln - 



|v/P ,„ / |v/|^ 



|v/lli^ 



GCT n 7 



f\\LUH-)\\y fhll + ln 



< ^||/|l^.^(^.)|lV/||Jl + ln 



Iv./H 
IvJil 

x,xJ II2 



1/2 



'csW^ixn' 



iv./ll 



dx 



where we have used the Jensen inequality - precisely, for a function g of integral 1, J gin 17 < hi(J g^) - in the 
last but one line, and the Sobolev imbedding from H^ into L* with constant Cg in the last one. The constant C 
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may change from line to line. Using this and Jensen inequality in the previous equation, we get 

\*WfMWf\m{u)dxdu< / \\MW fhW"^ fhmiu) du , 



<C'cT||/|lL2^(if2)||V/||2,.j||V/||2,A 1+ / - In - rfw 

Remark that the fraction inside the logarithm is always greater than 1 so that the square root is well defined. 
In order to get a bound on IJVj^/llj ™ , we use that 

I1V,/|1^ = 2 j h{u)^ = 2 j h{u)AiyJ^^du 

Va(|V,/n . — dxdu < 2||V|/||2||V,/||2 , (3.15) 

\u\ 

so that 

l|V./!l^_^ = !|V,/||2^„ + ||V,/||^ < 2 (||V,/||2,A + ||VaV,/||2,™) ||V,/||2,rn 

Therefore, 

1 



^dt\\VJh,n. + ^||V,,aV,/||^,^ < (41. + 2/3)\\^.f\\iuO + 

CcT|lVl^,/|l2,™|lV/|||,,^„ (||V./|||.^ + llVaV./|l|^) ( 1 + In 



c IIV^ /"IP 

2 



(3.16) eq:gevol3 



IIV./I12,™ 
We next use the inequality 1 + ln(a;) < — , valid for any e € (0, 1), a; > 0, and get 

dtWVJh^^ + i^WV^mVJWI^ < [Aiy + 2/?)|| V,/!^ „„ + . . . 
C(6)||V,,aV,/||i+^||V,/||}T^ (||V./|||,^ + ||V,,sV,/|||^ 

with C(e) = ^^.^cj. Now, with the temporary notations a — ||Va;/||2,m and b = \\yx,u'^xf\\2,m, what we need 
is to eliminate all the b in the right hand side, with the help of the b of the left hand side. Precisely, in the right 
hand side, we have the two terms 

b^+^a^-" , and b^+^J'" . 

We will use the Young inequalities xy < - — h — , where - + - — 1. For the first term, with p = -^ and then 
q = Y3J , and for the second with p — g— 2^ and then q — j^2i ■ ^^ obtain the following bounds 

7i4^f9-f .-! + £, 9 1^£ 4+ -2^ ,,3+f 3_, ^3 + 2en 1 — 2e c 1 st 
V+^a^ "■ < b^ H a^+i-' , and b^+^a^ " < b^ H a^+i-^. ^ 

valid for e < |. Taking into account the two constants ly and C{e) we get for e < ^ 

ldt\\\7j\\l^ + ^\\\/x,u^x.f\\l,n < (4^ + 2/3)||V,/||2 „ + . . . 

2 1 + e - 4-t- ^^ 4 3 + 2 e - fi-l ^^ 

C{e)^^-^\\yxfh^^-'+C{e)T^^-T^\\\7J\\^'^^-^^ . 
where C(e) has only change from a numerical constant. With the notation h — || V2;/||2 „, it gives 
-dth < {^v + 2(3)h + C(e) Ajz-^/i^+T^ + C(e)T^j/-T^/i^+T4fc , 

13 



That is a differential inequality with a growth faster than linear, that give a solution that may explode in a 
finite time. The time of explosion r* may be bounded below by something depending only on v, e, T and 
/i(0) = II V2:/i||2,m- Since e may be choosen arbitrairy between and 2 j t* depends only on i/,/3, T and the initial 
value /in- , ,., 

. . . . ffq • crPV0l3 

Using that in the inequality (p.iD|l . and turning back to the original u variable, we obtain the existence of a 
constant Cz{T,v, ||Va;/i||m) such that 

l|v./J|„, + ^y ||(v„a„)v,/|||,„di< C3(r,z.,||v,/,||2,,n) 

A bound by belov^r for r*. In order to simplify the next two paragraphs, we assume that jS <v. A careful 
analysis of the term in the right hand side of the previous equation shows that if e is choosen small enough so 

that V < C(e), then for /i < ft, := I ^^ 1 the dominant term is 4i//i and ioi h > h the dominant term is /i^+i-2« . 

Then, if h{0) > h, the explosion time is given by the equation 

dth < C(e)Tii/-T^/i3+T^ . 
For that equation, we get that the explosion time is larger than 

In the case h{0) < h, then for the early time, the equation may be rewritten 

dth < I2vh, 

till h{Q) = h. It take a time greater than T^ = j^ In ( -tot j . And after that time, the explosion time is given by 
the later calculation, and due to simplification it comes — . Finally, we get an explosion time 

^ . f h \ C 
■In 



12iy \h{0)J V 

Best choice for e. It is quite difficult to optimize that quantity in e. But, as the condition on e are < e < ^ 
and V <C{() = C^^^^e can choose 

e = mm -,C — 
\8 V 

With that choice, we get 

U(l + ln(-g)) else 



13 , — 



h) liv> 8CcT, T* 



jCi^^hg' iffto>l 



c 



^(1-ln/io) else 



It is then clear that the value of r* depends only on v, ho and the temperature T. 

The case of physical interest is the first one. Since in core of tokamaks, we have a large temperature T which 
implies a constant ct large, and a small colisionnality, in other words a small v. 
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: existence 



lem:approx 



4 Existence of solutions. 



existence 

-T' 

ent 



The proof will use the following notation and a preliminary 
ll'^llri'^ nlu.^' " ■ ----- 



^T] lead to the definition of 



In this section, we prove the existence theorem 

lemma. A priori estimates of the Lemma IbXl on'the solution (/, $) to| 

the set K of functions / such that 

||/(t)||2,m< VM, a.a.teO,T, 

where 

p2/3t _ -1 

M^WMlm + i'^^ + l^f-^WMh.u 

For each n > 0, we also introduce an approximation of the potential ^^ defined for any / G L'^ by 

I'.rt. t 



$„(i,x):= Y^ e* 

|fc|<n;fc^O 



1 - H{k) 



27rJ°/„(t, k,w)wdw — 1 



(4.1) 



lU} with 



Lemma 4.1. For any n G N* and any T > 0, there is a unique /„ in A'n L^(0, T; _ff„(f2)) solution to 

the potential $ replaced by $„ = $«(/«) and initial condition fi. That solution satisfy all the a priori estimate 

of the previous section. 

Proof ofpLcmma fa. lI Lct 6 be the map defined on K by S{f) — g, where g is the solution in KnL^{0, T\ H^{n,)) 
to llVip ' witn the potential $«(/) and initial condition /. The existence and uniqueness of S{F) follows from 
[[8j° Thnf T.T p 257, since V$„ is bounded in L°°{0,T;H^{T^)) by c„M for some constant c„. Then S maps 
K into K. Moreover, 5* is a contraction in L°°(0,T; L^(r2)) for T smallenoxigh. Indeed, let gi = 'S'(/i) (resp. 
52 — S{f2))- By estimates very similar to the one performed in Lemma IB .51 it noms that 

Vt > 0, ||(<i>„(/l) - $„(/2))(t, OIU-lT^) < C-„|l(/l - /2)(i, •)l|2,™, 

for some constant Cn- Substracting the equation satisfied by 52 from the equation satisfied by gi and integrating 
over il leads to 



g2(2^^^ 

2 di 
< 

< 
< 



^-2(2.+p)ty^ 



v\\{Vx:du){gi 
v\\{V:c,du){gi 

■ 2 



.92III 



52)112, m 



52; 



2 
2,m 



4zy2 



/2)ll2,mll52|l2,r^ 



< 



.g2V^(J°($„(/2) - $„(/i))) • V(.gi - 32) m{u)dxdu 

' /2)||2,m||V2:(5l - 52) || 2,m ||52 || 2,m 
(/l 



-c„||(/i 



4i^- 



/2)|l2,m 



And so 



151 



.g2||L~(0,T;L^) < cTe2(''+«^||/i - /2||l-(0,T;L?J. 

Hence there is a unique fixed point of the map S on [0, Ti] for Ti small enough. The bounds used for defining Ti 
being independent of Ti , a unique solution of the problem can be determined globally in time by iteration. The 
fact that this unique solution satisfy the a-priori estimates of the next section is clear since these estimates only 
depends on the bound satisfied by $ and not on its precise form. ■ 



existence 

Proof of Theorem 111. 1 

The sequence (/„) is compact in -^foc.«((Oj'^) ^ ^)- Indeed, it is bounded in 

L°°{Q,T\Ll{nu)) n L^{0,T;Hl{nu)) for any bounded subset U > (Recall of flu = T^ x M+). It follows 

IQ 

from the interpolation theory that (/„) is bounded in Lu ((0,T) x T^ x U). Together with the boundedness of 
C^x'^Mn)) in Lfnr „((0, T) X n), this implies that (^) is bounded in W,"^'f ((0, T) x flu)- By the Aubin lemma 

ir.i ons-Hagenes 'oc.ii at ioc,u 

re , It holds that (/„) is compact in L^{{0, T) x flu), so converges up to a subsequence to some functionfin^^^ 
It remains to pass to the limit when n — > +oo in the weak formulation satisfied by /„. A weak form of (|ll.l(l - (|ll.4|l 
is that for every smooth test function a with compact support in [0,T[xJ7, 

/ fi{x,u)a{0,x,u)udxdu 

nv^ xfn • ^ xCt + dufnduCt + f^u fnduCt ) dxduds . (4.2) eq:weakf on 



fn{-^ + V^(J°$„(/„)) • Vxa]u dxduds 
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J efljjreakf orm 

The passage to the Umit in fi.z\) wlicn n — > +00 can be performed if 

hm / f ufnV^{Jl{^n{fn))-'^xadxduds^ f f ufV ^ {Jlmf)) ■ V ^adxduds. 
"■->-°°Jo J Jq J 

This holds since (/„) (resp. {'^x{Ju{^n){fn))) strongly (resp. weakly) converges to / (resp. Va:(J°($(/))) in 
-^Lc uii^'"^) ^ ^)- ^^"^ since the /„ satisfy all the a priori bound, the limit / also satisfies them. H 



5 Short time uniqueness and stability of the solution. 

a^: uniqueness 



In this section we prove the shot time uniqueness and stability theorem 111. 1^ 

J lljjn.: uniqueness JaQ.igyroFP2D 

Proof of Theoremlll.:^! Ueiiotc by /i (resp. f2) a solution to lll.il for the field $1 (resp. $2), by Sf = fi — J2 
and by 8^ — J°($i ~ $2)- Multiplying the equation satisfied by (1 + u^)6f by 6f and integrating w.r.t. {x,u) 
with the weight u leads to 

Ij^W^fWlm < -'^IKV., a„)<5/||L„ + (4^^ + mWSfWlu + J SfV^iS'^) ■ V./2 m{u) dxdu 

To estimate ||(5/V((5$)||2.um, apply the inequality 

ab < e'' - b + blnb , a, 6 > 0, 

to (a, ^) = ( I giiy^g^ii ) , ( jsf\r I I fo'^ every nonnegative u and apply the Triidinger inequality (See \^^r~ 



I 



^\<^\\^pp\\2j dz<2. 



[lain.: rhor j fiffim: phi^ j 



Therefore, using ||V^(5$||2 < ■^^|lV5/||2,m (Lemmas Ilr5.5l and ||J.:^I and lb.4|) and the Jensen inequahty 



,: boundL 



(5.1) 



eq:exp-int 



WSf^Asn 



c\\v'sni\\sf\\ 



< C\\V'Smi\\5f\\Ul + 



6f y / \wj<^\ 

\SfhJ Ul|V2<5ci>||2 



dx 



■In 



\\sfh wm 



dx 



< 



ca 



|V.<5/||LJ|^/||^ l + ln 



\5fr 



U ■ ■- -V Vll'^/ll2 

Integrating in u with the weight um, it holds using again Jensen inequality that 



\6fyismlum < '2c4\\y.6f\\ij\sf\\i^J- 



\Sf\[ 



l^/lli,. 



< C4\\VJf\\lJ\Sf\\l^{l + ln 



\sn' 



1 + ln 

2 



ml 

1-5/111 



1-5/111,™ 



m{u) du 



Consequently, 



ld_ 
2di 



\\Sf\\l„,<CcT\\VJfh,n\\5fh: 



^ 



l + ln 



\sfr 



LUL") 



\sf\\L 



1V/2|12,^ 



^ll(V.,a„)<5/||2 +(4f. + 2/3)||<5/|| 



2,m > 
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App:A 



and finally using the inequality ||/(t)||i2 (^4) < e^*^ ^'^^'II/iIIl^ (^4) from lemma iB^T 

< -^¥f\\LM-j^^ jl|V/2||l:.+(4^ + 2/3)|l<5/|l^,„, 

where 

r = e^ (||/i,j||l^(l4) + II/2,j|Il?„(li)) ■ 
Defining s{t) = ;3j|l<5/|ll,„e-2W+2-)*, we get 

m<^\\^h\\Us(t)\n^ 
Av ' " s[t) 

It follows from the Osgood lemma that 

s(t) < s{Oy~"^'^ (5.2) |eq:osg 

with H(t) = -£^ L ||V/2(s)||2 ™. ds. We will show below that H is well defined on [0,t*], the time defined in 



ILeel: gxgaixl " 

Lemma b-Yl It implies 

||^/(i)||2,™<e(^+^^)^^-^"""||<5/(i)||-"'" 



\2.m 



and from that inequality we get the short time uniqueness and stability. Remark that the previous calculation 
do not use V^^/i and this is why we do not need an assumption on this quantity in the stability result. 



It remains to prove that H is bounded on [0,t*]. In fact, 



since (1 + u^)^ < 22 (1 + u^). Moreover, J^ IIV/2II2 „o dt may be bounded using Lemma h.(\ and the inequality 
(|l5.It)|l (we emphasize that the form used here is slighty different because we are in here in the ID setting for the 
variable u) 



W^fmiuodt < 2/ ||Vl„-/(t)||2,„||V./(t)||2,„ 

Jo 

< sup||V,/(t)||2,„V^( / \\ylsfit)h,udt 

t<T' \J0 J 



< C* 

V 



,: u-moment 



And /J" IIV/2II2 „7 dt is bounded by Lemma iB.a] with n = 3. 



A Appendix: Some controls on the Bessel function of zero-th order. 

The first Bessel function j" is much used in this paper. Indeed, in Fourier space, Ji that appears in the definition 
of the gyroaverage of the electric field, is the multiplication by J". Some properties of the function j" are given 
in |iuj. in this appendix, some bounds on J^ and its derivative are proven. 
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lem:boundJ Lemma A.l. J*^ satisfies the following estimates for all k C^M. 

.) |J°(.)|<min(l,^) 



m) |(jO)'(fc)|<mm 1,W^ 



2 

Trfc 



^) |(J")'(fc)|<(l + fc')"" 



J leiiLiboundJ 

Proof of Lemma IIA.ll 



First Inequality : The bound 17*^(^)1 < 1 is clear from the definition of J°, 

jO(^k)^^ f e''"'°''Ue = - f COs{k cos 9) d9. (A.3) |eq:JObis 

27r Jo Ti" ./o ' ' 

The bound by {^/2k)^2 is obtained as follows. J° is solution of the ordinary differential equation 

fc2(j")" + fc(j")' + A:V° = , J°(0) = 1 , (J°)'(0) = . (A.4) 

The new unknown u = ^/kJ^ is solution to 

There are no exact initial conditions for u. However, 

u{k) ^^[l + 0{k^)], u'ik) ^ J_[1 + 0(A;2)]. 
The second equation admits the fc-dependent energy, 

2 ' 2 V ' ik 



H{k) - H{k, u,u')^— + —[l + ^^ 



that satisfies 



H{k) - H{ko) =- I 

Jk. 



ko 4^3 



It follows from the behaviour of u near that 



H{k) ^ ± + 0{k). 
k->o+ 4fc 

Moreover the series expansion of j" near k — 0, 

^_k^ 



j°(fc) = E(-i) 



and its alternate character if fc < 2 imply that u^{k) > k — ^ (valid for k < V2)- Using, the inequality and the 
behavior of H near 0, we get if < fco < /c < \/2 

1 „,. . /"V 1 1 






since the first line is satisfied for any fcg > 0. Therefore, 



u^k) < fc^, k<V2. 



18 



A simple calculation shows that the function appearing in the right hand side is increasing in fc, so that 

u\k)<^, fce[0,V2]. 

/— J egiiik 

For k > v2, simply remark that H is decreasing and that from (|IA.b|) 



u^ 



(fc) < 2H{k) < 2H{V2) < -^ 

v2 



In any case we get v?{k) <2 2 which gives the desired inequality. 

Second inequality : It is a consequence of the first, for fc > 1. For fc < 1, it may be obtain from a comparison 
of the entire development of J^ and (1 + fc)^^/'' around the origin. We get 

1.2 1,4 l2 c?,4 ^r.h^ 

jO(fc) <l_L + L<l_L + £L_i2^<(l + fc2)-l/4 

^ ' - 4 64 - 4 32 128 - ^ ' 

| pr|- . Tn 

Third inequality : Taking the derivative of J in the definition (pi.a|) . 

(jO)'(fc) = A / coseie^'^^'^de^-i / cos6'sin(fccos6')d6', 
from which it is clear that |(o/^*')'(fc)| < 1 for all fc. Next we transform the previous integral in 

(J°)'(fc) = -^['^^^Mda, 






V< 



where (/ii)i<i<j are the points where sin(fc0) vanishes and 1, 

TT , 27r , (7 — l)7r 

/i0 = 0<ft.i = -</i2 = -;-<■■•< /ii-l = , < hj = 1. 

k k fc 

The previous sum has alterned sign, the larger terms occuring for large i. Its terms are with increasing absolute 
values, except for the last one which is incomplete and may be smaller than the next to last term. However, 

i 

-Sl < So - Sl < '^{-lySj < So - Sl + S2 < So, 
i=0 

so that 

|(J°)'(fc)l < max(so,si)< ^ ^ "'^^ 



I" Jl-7r/fc VT 



1 /27r 7r2 / 2 



This ends the proof of the third inequality. 

The proof of iv) is similar to the proof of m), since J ^ < 2" 
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